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Abstract 

We study flux penetration in a disordered type II superconductor by simulations of 
interacting vortices, using a Monte Carlo method for vortex nucleation. Our results 
show that a detailed description of the nucleation process yields a correction to the 
scaling laws usually associated with flux front invasion. We propose a simple model 
to account for these corrections. 
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1 Introduction 



The magnetization of type II is traditionally described in terms of the Bean 
model (|l|): magnetic flux enters into the sample from the boundaries, form- 
ing a flux gradient that is pinned by disorder. At the microscopic level, the 
process takes place through the nucleation of vortex lines carrying each a flux 
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quantum The Bean model provides a phenomenological picture of average 
magnetization properties, such as hysteresis and thermal relaxation (j^, but it 
is inadequate to describe fluctuations, which turn out to be quite important. 
As it has been shown experiment allyusing magnetoptical methods, flux fronts 
are typically rough or even fractal 13; S 0)- 

Recently, we have shown that the flux penetration in a disordered supercon- 
ductor can be described by a disordered non-linear diffusion equation 
The equation can be obtained performing a coarse-graining of the microscopic 
equation of motion of the vortices. In the absence of pinning, it reduces to 
the model introduced in Ref. (Q). This model has been solved analytically to 
provide expressions for the dynamics of the front for different boundary con- 
ditions (9^ .1Q |) and the results are in perfect agreement with vortex dynamics 
simulations (jsj). When quenched disorder is included in the diffusion equa- 
tions, flux fronts are pinned in agreement with individual vortex simulations. 
Varying the parameters of the equation, we observe a crossover from flat to 
fractal flux fronts, consistent with experimental observations. The value of 
the fractal dimension suggests that the strong disorder limit is described by 
gradient percolation (0). 

Here we reconsider the influence of boundary conditions in the front dynamics. 
In Ref. (0) we have analyzed several types of boundary conditions determining 
the way vortices enter into the sample depending of the experimental setup 
one would like to model. For instance, a constant applied magnetic fleld can 
be approximated by a constant vortex density at the boundary of the sample 
13 nu \. This assumption simplifles the real vortex nucleation process and 
allows for straightforward numerical and analytical approaches (|7j; 0; [l3 0)- 
A more accurate description of vortex nucleation can be obtained combining 
the vortex dynamics simulations with the Monte Carlo method (0). Using 
this method we find that the simple widely employed approximations for the 
boundary conditions are only asymptotically true. At short times the details 
of the nucleation process affects the expected scaling behavior for the front 
dynamics. We are able to quantitatively estimate the corrections to scaling, 
using a simple front dynamics model in the spirit of Washburn approach to 
fluid imbibition (E^ 



2 Simulations 



We use here the model introduced in Ref. (jlll). which combines a typical vortex 
dynamics simulation scheme with a Monte Carlo method for vortex nucleation. 
In a very large sample with a constant magnetic fleld H oriented along the 
z axis, vortices are modeled as a set of interacting particles performing an 
overdamped motion in the xy plane. 
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The Gibbs potential associated to vortices of coordinates fi can be written 



as 



ij i i 



where in an infinite system Ei2{r) = [$Q/(87r^A^)]i<'o(M/'^) is the vortex- 
vortex interaction in the London-London theory, with fiux quantum $o and 
penetration length A. Here we consider a semi-infinite system, bounded by the 
y = line, and consequently we add to the vortex-vortex interaction a term 
accounting for the interaction between each vortex and the image of the others 



( 111). The term rj(xj), where Xi is the distance between the vortex i and the 
sample surface, represents the interaction between each vortex and its own 
image Finally the external magnetic field gives rise to a sort of chemical 
potential with = $o(l — exp(— x/A). In Ref. (fill ) the interaction energy 
with random pinning centers is also included to Eq. 1, while here we restrict 
our attention to a clean system. 

The vortices evolve according to an overdamped equation of motion Tdtfi = 
—dr-G, where F is a damping constant. The equation of motion is integrated 
numerically and after each integration step a zero temperature Monte Carlo 
step is performed: a new vortex is nucleated at a random position in a strip 
of length A close to the sample boundary if the Gibbs potential is reduced. In 
practice, we consider a system of size Ly = 20A and = lOOA, with periodic 
boundary conditions along the y direction. At the beginning of the simulation, 
we start with an empty lattice and vortices are nucleated close to the a; = 
boundary. Since we are only interested in the transient behavior the boundary 
condition at X — Lx is inessential. 

As the vortices are nucleated, they are pushed toward the interior of the sample 
giving rise to a density profile. The evolution of the profile is reported in Fig. 1. 
The results are in slight disagreement with the simplified boundary condition 
used in Refs. ^ E3 @)- ^ constant boundary vortex density. Fig. 1 clearly 
shows that the boundary density increases. In addition. Fig. 2 shows that the 
front position does not evolves as a power law, Xp ~ t^^^, as predicted by the 
theory JlliSli). 



3 A simple model for front dynamics 



In order to account for the behavior observed in numerical simulations, we 
consider a simple model for the front dynamics (@). The front is driven by the 
density gradient which can be can be estimated simply as Vp ~ Po/xp, where 
Xp is the front position and po = p(0, t) is the boundary density. Thus the 
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Fig. 1. The density profile at different times t obtained from numerical simulations. 
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Fig. 2. The position of the front at different times t obtained from numerical simu- 
lations. The dashed line is the t^/^ prediction 
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Fig. 3. A comparison between the front dynamics model and the numerical simula- 
tions. 

equation of motion of the front can be written as TdXp/dt = ap^/xp, where 
a is the interaction parameter computed in Ref. (j^. In order to solve this 
equation, we have to specify the dynamics of the boundary density, which to 
a first approximation can be written as rdpo/dt = Ph — Po where pn oc H is 
the asymptotic value of the vortex density, and r a characteristic time. 

Solving the two differential equations, we obtain 



Po = P/f(l-exp(-t/r)) (2) 
Xr>{t) = (^)'^' {t + r(l - exp(-t/r)))V2. (3) 

In order to compare this result with numerical simulations, we plot on the same 
graph the numerically calculated po and dXp/dt together with the theoretical 
prediction from Eq. 2. The numerical result can be well fit by the model with 
T ~ 10. 
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4 Conclusions 



In this paper we liave considered the effect of vortex nucleation on fiux front 
propagation in type II superconductors. We have studied the case of a constant 
applied magnetic field and compared the results of numerical simulations with 
previous approaches to the problem While asymptotically we recover 

previous results, a more accurate account of the vortex nucleation process leads 
to important corrections for the scaling laws. This fact should be considered 
in the interpretation of experimental results. 
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